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ABSTRACT 

The issue of developing effective and robust schemes to implement a class of the Ogden-type hyperelastic 
constitutive models is addressed. To this end, explicit forms for the corresponding material tangent 
stiffness tensors are developed, and these are valid for the entire deformation range; i.e., with both distinct 
as well as repeated principal-stretch values. Throughout the analysis the various implications of the 
underlying property of separability of the strain-energy functions are exploited, thus leading to compact 
final forms of the tensor expressions. In particular, this facilitated the treatment of complex cases of 
uncoupled volumetric/deviatoric formulations for incompressible materials. The forms derived are also 
amenable for use with symbolic-manipulation packages for systematic code generation. 

1. INTRODUCTION 

To a large extent, constitutive models of the so-called generalized Rivlin-Mooney type [1,2]; i.e., 
with the (stored) strain energy density written as a polynomial function in terms of the deformation 
invariants , have dominated the phenomenological theory of isotropic hyperelasticity [1-6], as well as the 
related computational literature on finite-strain elasticity [7-9], over the years. Recently, alternative 
representations in terms of the principal stretches have become increasingly popular in nonlinear finite 
element analyses [6,8,13]. In particular, several forms in this category were given by Ogden [5], Peng 
[11], and Valanis and Landel [12]. From both the mathematical and physical standpoints, these latter 
models exhibit a number of attractive features. For instance, in the Valanis-Landel hypothesis [12], also 
called the separability assumption, the underlying their mathematical structure (i.e., a strain energy 
depending on separate symmetric functions in each of the principal stretches) was found to be essential 
for effective and accurate characterization of the behavior of incompressible, and slightly-compressible, 
rubbers at (relatively) high strain levels [5,6]. 


However, from the viewpoint of numerical implementation, the use of these models presents a 
number of unique and difficult problems, which do not arise in alternative representations in terms of the 
strain invariants. Basically, the main difficulty here stems from the fact that, in addition to their being 
“reasonably complicated" functions of the strain components, the orthonormal basis of eigenvectors 
associated with the principal stretches (eigenvalues) are not uniquely defined for the whole deformation 
spectrum; i.e., nonuniqueness for the case of two or three repeated eigenvalues [14,16-18]. This will be 
a major source of difficulty in obtaining the (directional) derivative of the tensor-valued stress function 
(in terms of deformation tensor) needed to obtain the elastic-moduli matrix (i.e., material tangent stiffness). 
These latter arrays are essential ingredients in any nonlinear incremental finite element solution; e.g. using 
either the total Lagrangian , or Eulerian or “updated" Lagrange approaches [8,20]. General solutions of 
this type will typically involve very complex, nonproportional deformation paths during which the number 
of uniquely-defined eigenvectors will change, and the necessity of obtaining explicit forms for the material 
stiffness in such singular or limiting cases becomes obvious. 

The type of problem mentioned above is not new; it has a long standing in various applications in 
continuum mechanics [21-28]. An important case to the point concerns the use of Hencky-type, or 
logarithmic strain measures in various contexts [25,27]; e.g. in conjunction with derivations of alternative 
forms for its (material) time derivative [22,24], the associated conjugate stresses [25], the so-called Hill’s 
principal-axes techniques in the numerical treatment in [26], etc. However, to the present authors’ 
knowledge, the topic has not been previously addressed in the context of deriving material-moduli tensors, 
for the general three-dimensional analysis. For suitable particular forms in the case of two repeated 
eigenvalues, but restricted to the two-dimensional problem, we refer to [8,13] for simple derivations. 

Addressing the above problem constitutes the main objective in the present paper. In particular, we 
give details leading to the derivation of the proper explicit forms for the stress-deformation relationships, 
and the associated material elasticity arrays, in terms of general components of the deformation tensor for 
the entire regime of distinct as well as repeated; i.e., double or triple coalescence, eigenvalues of the strain 
tensor. In this, two specific forms of the Ogden-type strain-energy functions are considered; these are of 
sufficient generality to include most (if not all) of the models currently in common use. The first form is 
for the so-called general, unconstrained, materials using the traditional definition of principal stretches. 
The second is suitable for treating incompressible, and slightly compressible, cases and it makes use of 
a modified set of principal stretches. These are associated with the concept of deformation decomposition 
into pure dilation-isochoric distortion response components, which has been utilized extensively in recent 
literature [5,8,9,13], including works on finite elastoplasticity [9,10]. For both of these models the 
implications of the underlying separable property are exploited throughout the derivation, thus providing 
a convenient means for performing the limiting process. For instance, this eliminates the need for a 
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posteriori search and identification of various “singular” expressions that may (appropriately) combine to 
give nonzero limits in scalar/tensor-valued composite terms [8,13]. 

2. PRELIMINARIES 

For convenience, we summarize in this section some fundamentals and results pertinent to the large- 
deformation analysis. These will be referred to repeatedly in subsequent discussions and derivations. Here, 
and for the remainder of the paper, both indicial and the “counterpart" matrix notations will be used 
interchangeably; repeated indices imply the summation convention. Only rectangular Cartesian coordinate 
systems are employed, and we use lower-case subscripts for all tensor indices. For the most part, this is 
sufficient for our purpose here, and there is no need to consider lowering or raising of indices. Note, 
however, that the only exception to this concerns the definitions of “pull-back” and “push-forward” of 
tensors given below. 

2.1 Basic-Kinematics 

In line with common usage [e.g. 14-16], the set-up used to describe the nonlinear kinematics of a 
continuum is as follows. Let P be a material point with position vector X in the undeformed (initial) 
configuration “B”. The map of P in the (current) deformed configuration “b” is denoted by p whose 
location is given by x. The deformation map is defined by (J): B — > b, and 

jt=(J)(*,r) ; F = a 4>/aAT (or F ij = dfy/dXj) (2.1) 

where F is the (two-point) deformation gradient tensor for the current configuration “b” (at time t). 
Because of the one-to-one nature of the deformation map, F is nonsingular and it can thus be expressed 
as follows, using the (right) polar-decomposition theorem, 

F = RV ; or F,j.R it U kJ (2.2) 

where U is the symmetric, positive definite, pure stretch tensor , and R is a (proper-orthogonal) rigid- 
rotation tensor, such that 

R R t = / ; C = U 2 = F t F ; J = det U = det F ( 2 - 3 ) 

where I is the second-order identity tensor, the superscript “T” indicates a transpose, and C is known as 
the (right) Cauchy-Green deformation tensor and J represents the change in volume. 

Through eqs. (2.1)-(2.3), several other strain measures can be defined. For example, we recall here 
the following familiar expressions for Green-Lagrange strain, E, and Almansi strain, e, respectively, 
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E = - (C-I) 

2 


(2.4) 


e = .!.(/ - F~ T F~ l ) (2.5) 

where ( )'* indicates an inverse. These are material and spatial strain fields [Sec. 2.2], respectively. 

As an alternative to Eq. (2.2), we may employ the following decomposition for F: 

F = J m F ; C = J~ 2/3 C = F t F ( 2 - 6 > 

1 /O A /V 

into its two constituents; i.e., pure dilatation , J I, and isochoric distorsional deformation, F (det F = 

C 

1). Correspondingly, we refer to the Lagrangian tensor as the modified , volume-preserving, deformation 

tensor. This provides for a proper treatment of uncoupled volumetric/deviatoric response of, for example, 
incompressible rubbers [5,8,13]. 

Considering the second of Eqs. (2.6), the (directional) derivative formula [14,15] for C leads to: 


DC ■ A = J 




for any variation A in tensor C, and we thus have 


(2.7a) 


dC/dC = J 


- 2/3 


/ ( 4 ) -Ic® 

3 



(2.7b) 


where symbols : and <8> denote scalar multiplication (i.e., trace operation) and vector product of tensors, 
respectively; and the components of the fourth-order, symmetric, unit tensor I*- 4 * is given by 


'm - W 


( 2 - 8 ) 


2.2 Material (Convected) and Spatial Description 

The geometric notions of “pull-back" and “push-forward" operations for tensor and vector fields 
on manifolds [15] play an important role in developing “counterpart" expressions for the stress- 
deformation and material-moduli (elasticity) tensors in both material (Lagrangian) and spatial (Eulerian) 
settings. They provide the appropriate forms in finite element formulations based on total and updated 
Lagrange methods, respectively [8,20]. 
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Proper definitions of these operations require identification of the precise placement of tensor 
indices (i.e., “covariant" or “contravariant" components), even in a rectangular Cartesian setting. Therefore, 
considering the case of “contravariant” tensor (vector) fields that are (independently) defined on various 
configurations, and letting T be a second-order material tensor field defined on “B”, its push-forward (J^T 
is a spatial tensor, T', of the same contravariant type on “b", or <j>(B), which is given by 

t' =<t> = ftf t ( 2 - 9 ) 

Alternatively, letting t be a “contravariant" two-tensor field on “b", its pull-back, t', on “B" is defined as 

t * (j)*< = F 1 t F t (2 ‘ 10) 

As an application of this latter operation, we recall here that the (symmetric) second Piola-Kirchhoff 
stress S is actually the pull-back of the Kirchhoff stress r = J a (where a is the Cauchy or true stress 
tensor); i.e., 

5 = = F~ l T r r (2 - n) 

Similarly, considering the (covariant) components of the material and spatial strain tensors E and e, we 
have 


e = 4>*£ = F~ t EF' 1 


( 2 . 12 ) 


Finally, we note that the above notions are also useful for deriving different forms for objective 
rates for spatial stress and strain fields [14-16]. For example, we cite here the two most commonly used 
rates of this type; i.e., the Truesdell rate of Cauchy stress, a, and the rate of deformation tensor, d(t = 
time, and an overdot signifies a material time derivative): 




— fsf t 

J 


(2.13) 


d 



= F~ T E F~ x 


(2.14) 


In view of the above equations, we actually proceed for the remainder of the paper with derivations 
in terms of the material-type tensors S and E and the associated rates; i.e., for the stress-deformation 
relations and material-moduli (fourth-order) tensors D (recall Eq. 2.4): 
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S = S(C) ; S = £>(C)£ 


(2.15) 


where ( ) indicates “a function of". The “counterpart" spatial forms follow directly by the push- 
forward/pull-back according to Eqs. (2.1 1) to (2.14). In particular, for the spatial tensor D one has 

Jt,.TKe)J ; JV - F im f>, C ) V„ (2.16) 


2.3 Spectral Representation 

Denoting by X i and n- (i = 1,2,3) the principal values of C (i.e., square of principal stretches) and 
their associated eigenvectors, respectively, we can utilize the following representation for its components: 

/= 1 M 


or, equivalently, in terms of the (rank-one) 3x3 matrices N i? 

(i) for the case of distinct X^: 

C = £ \ i N i (2- 18a ) 

i=l 


(ii) for the case of double coalescence (Xj. * X s = \ = X): 

C = (\ r -\)N r +\I (2.18b) 

(iii) for the case of triple coalescence (Xj = X 2 = X 3 = X): 

c-xEjv,.-xz (2 - 18c) 

i=i 


In the above, the (r,s,t) is any cyclic permutation of (1,2,3), and tensors Nj are often referred to 
as the (orthogonal) eigenprojection operators (i.e., on the null space of C - Xj I); e.g., see [18,19] for 
further details. Note that, among the latter tensors, only those which are uniquely defined (i.e., 
corresponding to distinct eigenvalues) are employed in Eqs. (2.18). 
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Explicit expressions for Nj in terms of C can be obtained by suitable manipulations of Eqs. (2.18). 
This leads to the following forms corresponding to previous cases (i) and (ii): 


(0 


(ii) 


iV r = 


(X r -X s )(X s -X t ) 


(c-yxc-y)] 


N r = 


(X r -X) 


( C-XQ 


(2.19a) 

(2.19b) 


These will prove very useful in obtaining the derivatives of the stress function S = S(C) in Secs. 5 
and 6. 

Remark 2.1 . Alternatively, following standard procedures: e.g. as outlined in [17, p. 563], we can show 
that 


(K r -X s )(X r -\ t )N r = Cof (C-X r I) 

iL 

where Cof ( • ) denotes the adjoint matrix whose elements are the cofactors of ( • ); i.e., with its (i,j) 
component given by 

" 2 " ^mpi £ nqj ^ mn ^ r ® mn ) ^pq ^ r ^ pq ) 


where e is the rank-three alternating tensor. Together with the identity: 


this Cof tensor can be simplified to exactly the same form in the bracketed term on the right-hand side 
of Eq. (2.19). 


3. STRESS-DEFORMATION RELATIONS, THE SEPARABLE AND UNCOUPLED 
VOLUMETRIC /DE VI ATORIC FORMS FOR STRAIN ENERGY FUNCTION 

Restricting attention to the hyperelastic isotropic material case, the two forms of the strain-energy 
function, W, considered here are written as: 

w - W(X t ) = £ a tl (x“" + + X 3 ") (3-1) 

n = 1 


and 


W = W(\)+h{J) 


(3.2) 


where X are the principal values of C (see Eq. 2.6); h is a (convex) function of one variable J = (Aj A 2 
X 3 ) 1/2 , and the a n and a n are material constants. The functional dependence of W in Eq. (3.2) is assumed 
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to be exactly in the same separable form as in Eq. (3.1), but with the modified stretches X r now replacing 
X i (where = J" 2/3 \ from Eq. 2.6). 

Note that in most practical applications, the h(J) is actually taken as a penalty function imposing 
the incompressibility (or rather, the slight-compressibility) constraint [5 ,6, 8,1 1,12]. In this connection, we 
also note that an implicit assumption of strict incompressibility is often made [5,6] when using Eq. (3.1). 
However, such an assumption is not implied here, and for emphasis we refer to this material model as the 
unconstrained general case. Of course, adding a penalty-like term for the purpose of practical usage is the 
same as for Eq. (3.2). However, the separate treatment of Eq. (3.1) will provide the basis for subsequent 
derivations of the D tensor in the more complex situation of Eq. (3.2). 


The constitutive functions S = S(C) and D = D(C) can now be obtained 


or 


S 


2 — 

dC 




D 


A ^ 

dC dC 


dW d£ + h , dJ_ 

~dC~dC + dC 



(3.3) 


(3.4) 


where the primes denote differentiation of h with respect to its argument, J. Together with Eq. (2.7), the 
expression in Eq. (3.4) can be reduced to 


S 


= J 


-2/3 


DEV 


\ 

2 aw 

acj 


+ pJC~ l 


(3.5a) 


with 


DEV (B) = B -1(B:C) C' 1 


(3.5b) 


for any Lagrangian (stress-type) variation tensor B, and where p = h' is the (hydrostatic) pressure function. 
The DEV operator in Eq. (3.5b) allows for a concise writing of several expressions given later, physically 

it gives the deviatoric part of the material tensor A, see Sec 5.2. Its spatial counterpart: e.g. dev 
a = a - — tr (a)/, where tr( • ) indicates the trace of ( • ), is an already well-known operation. 
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4. MATERIAL-MODULI TENSOR. THE GENERAL UNRESTRAINED MODEL 


4.1 Continuous Representation of the Constitutive Stress Function 

We start by first considering the case of Eqs. (3.1) and (3.3). Making use of the fact that, for the 
present isotropic material, S is coaxial with C (he., identical eigenvectors for both), it admits a 
representation in terms of its principal values S i (i = 1,2,3), and directions N; (Eq. 2.19), in exactly the 
same form as in Eqs. (2.18). Explicitly, for the present separable W-form in Eq. (3.1), we have 


i.e. the same functional dependence S on only for each separate i = 1,2,3, where a comma-subscript 
with W indicates (partial) differentiation with respect to (w.r.t.) the following principal-stretch 
“coordinate”. Combining Eqs. (2.18, 2.19), and (4.1) leads to 



(4.1) 


(i) S = a C 2 + b C + c I 


(4.2a) 


(ii) S = a C + b I 


(4.2b) 


(iii) S = S(X) I 


(4.2c) 


where (i)-(iii) correspond to the three cases considered previously in Eq. (2.18), and 
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a = -£ m(X r -X s )S(k t ) 


(4.3a) 


r=l 



(4.3b) 


3 


* = "E nt(X r -X s )X,X s S(X t ) 


(4.3c) 


m 


(X { -X 2 ) (X 2 -X 3 ) (X 3 -X l ) 


(4.3d) 



(4.4a) 
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XS(X r )-X r S(X)]/(X r -X) 


(4.4b) 


b = - 


Recall that the (r,s,t) are still cyclic permutations of (1,2,3), and all notations follow those used previously 
in Eqs. (2.18). 

With the above expressions, one can easily proceed to investigate the continuity of the functional 
representation of S = S(C). In fact, this simply amounts to showing that the limits for Eq. (4.2a) exist, and 
that these are given by Eqs. (4.2b) and (4.2c) when the distinct eigenvalues (A^A^Af) approach (X r ,X,X) 
and (A, A, A), respectively. For example, considering the former case, we let A= (A r ,A 5 ,A r ) — > A = (A r ,A,A) 
along the line normal to the plane A s =A t , and there results, for the limits of a(A), b(A), etc. 


Lim a = ! [S (X r ) - S(X) - (X r - X)S'(X )] 

X - X (X r - X) 2 

Lim b = ! !(xj - X 2 ) S' ( X ) - 2X [S ( X r . ) - S(X ) 

X - X (X r - X) 2 1 


Lim c = 


1 


h (x r -xy 


X 2 S(X r ) + X r (X r -2X) S(X ) -X,X(X r -X) S'(A)] 


Lim C 2 (X) = 1 

X - X 


(X r ~X) L' 


x 2 r - X 2 ) c( X ) - (xxj - X 2 xJ I 


1 2 _ \ 


(4.5a) 

(4.5b) 

(4.5c) 

(4.5d) 


Recall that a prime indicates differentiation w.r.t. the shown (single) argument. Substituting Eqs. (4.5) into 
(4.2a), we obtain the expression in (4.2b). 

Similar arguments can be used for the case of triple coalescence, by finding the limits of S in Eq. 
(4.2a) when distinct X = (X r ,X s ,X t ) -* X = (X,X,X) along the line normal to the principal space diagonal, 
or hydrostatic axis [e.g., 17], X,. = X s = X t . The result is expression (4.2c). 


Remark 4.1 . Based on available results from strict mathematical analyses; i.e., the so-called Ball’s Lemma 
[28], the continuous differentiability of any tensor-valued function like S = S(C); e.g. its elastic moduli- 
tensor D of concern here (see Eq. 3.3), is guaranteed if we can simply show that the limits of D(X) with 
distinct X = (X r ,X s ,X t ) exist for the “degenerate" cases when X(X r ,X,X) or X ->(X,X,X) along the two lines 
that are normal to the plane X s = X t , and to the line X,. = X s = X t , following similar arguments as the above. 
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The presence of an individual \ -argument for each SCA^-function in various terms of D(C) will provide 
for considerable simplifications here. 

4.2 Explicit Expressions for D 

The general expression for tensor D can be obtained by applying the directional derivative formula 
to Eq. (4.2a), recalling the second of Eqs. (3.4), and making use of the formula dX JdC - N r (see 
Eq. 2.17) for the present case (i) of distinct Aj(i=l,2,3). Although tedious, it is straight-forward to show 
that 


where we define, for any two (second-order-symmetric-tensor) arguments G, H, the components of the 
fourth-order tensors, P and Q, as 




(4.6a) 


+ a 5 lQ(CJ) * Q(I,C) * } * 2a/ 4) 


P ijkf G ' H ) = G ik H jl + G il H jk 
Qijkf G-H) = G ik H jl + G il H jk + G jl H ik + G jk H il 


(4.6b) 


(4.6c) 


and r 4 -* is given in Eq. (2.8). We have also utilized the following definitions for scalar coefficients 


^ (i=l-6): 


3 3 3 


“r E 1r ; a 2 * E (.K-h^r i a 3 = E h V X r 


(4.6d) 



3 3 


«4 = E -K) 1 nr ; «5 = E + 


(4.6e) 
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H = E n r (W 2 + (V- / i)M r 


(4.60 


in which 
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/j Aj X^ X/j \ I 2 XjX2 ^2^3 XjXj 9 Ij - XjX^X^ 


(4.6g) 




5'(X r ) + (X, - X r ) (M r + M,) + (\ ~ \) (M r + M 5 )] / [(K - x/(X r - X,) 2 ] ( 4 - 6h ) 


= 5(X r ) / [(X r -X s )(X r -X,)] (4.6i) 

where the principal values S(\) and their derivatives w.r.t. X r , S'(X r ), are found from Eqs. (4.1), and the 
(r,s,t) is again any cyclic permutation of (1,2,3). 

Considering now the two degenerate cases (ii) and (iii) of double and triple coalescence of \ 
values, respectively (Eq. 2. 1 8), we proceed in exactly the same manner as in Sec. 4. 1 to find the associated 
limits of the expression in Eqs. (4.6). Omitting all the details and considerable simplifications involved, 


we summarize the results as follows: 

(ii) D = b x P{C,C) * b 2 [ Q(C,I) * Q(IC)] + 2 V (4) 

(iii) D = 2 S'(X)/ 4) 


where 


b x = 


a 


(X r - X) [s'(X r ) + S'(X) ] - 2 [s(X r ) - S(X)| 


(4.7) 

(4.8) 


(4.9a) 


b 2 = a{(X - X r ) [X 5'(X r ) + X^S'(X)] + (X r + X) [s(X r ) - 5(X)] } 


(4.9b) 


b 3 = a 


X r -X) 


X^ S'{\) + X 2 S'(X r ) 


•2X r X 


[5(X r ) ~ S(k) 


(4.9c) 


a = l/(X r - X) 3 ( 4 - 9d ) 

Remark 4.2 . (i) All formulas in Eqs. (4.6) - (4.8) exhibit proper symmetries in the D arrays; i.e., with 
symmetric strain rates £ (and C), and stress rate § tensors (see Eq. 2.15), = 1/2 ( dS- / dCy 

+ dS- / dCfc ). (ii) In view of the existence of the collective limits in Eqs. (4.7) and (4.8), the continuous 
differentiability of D follows immediately from results in [28] (although this is not true for the constituent 
arguments, and and Nj of C). (iii) Restricting the dimension of the underlying space of all tensors in 
Eqs. (4.6) - (4.8) to two, it can be shown that the D arrays corresponding to the general distinct- 
eigenvalues case (e.g., = X r * X 2 = X) and the repeated-two-eigenvalues case (i.e., X A =X 2 = X) are 

actually those in Eqs. (4.7) and (4.8), respectively; no separate treatment is therefore necessary for this 
two-dimensional problem. 

3 1/2 

Remark 43 . Being also coaxial with C, the logarithmic strain, fin U = (fircA r )N r , for distinct X r 
can be treated similar to the present case; i.e., its (material) time derivative r fikes the forms of Eqs. (4.6) - 
(4.8), except now its principal values (in parenthesis in the above summation) replace S(X r ). This may 
prove useful in, for example, large-strain elastoplastic finite element applications as in [e.g. 27]. We also 
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note that the present formulas for the spectral representations of “coaxial” tensor-valued functions (e.g. 
S here), and their derivatives, in terms of the argument tensor (e.g. C here), will prove very useful in 
applications involving continuum-based damage constitutive models and the experimental characterization 
of material coefficients, e.g. see [30]. 


5. MATERIAL-MODULI TENSOR. THE UNCOUPLED VOLUMETRIC /DEVIATORIC MODEL 

The analysis in the present case proceeds parallel to that given in the previous section. We 
therefore draw extensively on the results of Sec. 4 and omit many details. 

5.1 Continuous Representation of the Constitutive Stress Function 

In view of the coaxiality of tensor C, S, and C, and making use of Eq. (3.2) to explicitly evaluate 
the term DEV in Eq. (3.5), one can conveniently write the stress function in the present case as: 

S - J _2/3 3 +/(X,) C~ l +pJC~ l ( 5 - 1 a ) 

with 

8 = 2 dW/dC ; S r = 2W <r = S(k r ) ( 5 - lb ) 

M)--| E *,■*, < 51c > 

J r=l 


In this, the S r (r=l,2,3) are the principal values of S,, and the comma subscript in W defines 
differentiation w.r.t. modified parameters X. We make the important observation that, with Eq. (3.2), theS r 
displays exactly the same symmetric-single argument dependence as in Eq. (4.1); i.e., S admits the same 
limiting forms as in Eqs. (4.2) in terms of C. . 

With the neat separation of its constituents, no further elaboration is actually needed regarding the 
representation of S in Eq. (5.1a). The last two terms on the right-hand side are always well-behaved for 
any X;, whereas the continuity of S for distinct as well as repeated X; follows immediately from the 
analysis in Sec. 4 (simply replace Xj by X i; , C by C,, and S by S in pertinent paragraphs). Note that this 
will also be the case for constituents of the D tensor given in the sequel. 

5.2 Explicit Expressions for D 

A careful examination of the derivatives of various terms in Eq. (5.1a) w.r.t. components C reveals 
that, for all cases (i.e., distinct as well as repeated eigenvalues), a single general expression can be written 
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for D; some of the terms will, however, need to be “branched out" for their individual explicit evaluation. 
That is 

D = J~ 2r3 D§ + Y + (p"J+p'J 2 )C~ l ® C l - 1$ 0 C~ l + 2C' { ® DEV (M) (5.2a) 


in which 


Y = 




(5.2b) 


where the primes in p' and p" indicate first and second derivatives, respectively, w.r.t. J, the tensor 
operator? is defined as in Eq. (4.6b), and the newly-introduced tensors DS (fourth-order) and M (second- 
order) are defined by 

D§ = 2 d§/dC = J ~ m A (5.2c) 

3 


M = df(Xj)/dC = M(C) 


(5.2d) 


The above decompositions in formulas (5.2c) and (5.2d) are made to facilitate the use of previous 
expressions; namely, by its definition the (fourth-order) tensor A = A(C) is obtained in exactly the same 
manner as was for D(C) in case of the unconstrained model of Sec. 4. Thus, using the following 
replacements (assignment “operators”): 

C-C ; k^X, ; S(X)<-S(X) ( 5 - 3a ) 


and 


X r ) for i = 1 to 6 (5.3b) 

b i (X r )+-b i (X r ) for i = 1 to 3 ( 53c ) 


we can now write 

(i) A = Eq. (4.6) ; (ii) A = Eq. (4.7) ; (iii) A = Eq. (4.8) 


(5.4) 
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Similarly, we can use the above assignments in finding the scalar product of A and C needed in 
the last term of Eq. (5.2c); see the scalar operator C‘* : A of Eq. (2.7a). There results, for the (second-order) 
tensor B, and upon use of the Cayley-Hamilton theorem [14,15], 

g(m) _ q + q 3 ' 2 I ; for m = (i) or (ii) 

B m = S(i)C (5,6) 

where for Case (i) : 


< h 


h(h~ 2 h) + h 


+2(/ 1 -/ 2 )a 2 + h( 2a 3 +a 4) +2a 5 


;2 t 2 


q 2 ~ 0\h + ^2 ~h^2) a \ + 2 ( 3 hh^ a 2 h( 2a 3 + a 4^ +a 6 


<jf 3 - / 3 (/j ~ I 2 )a i + 2/|/ 3 « 2 + / 3 (2a 3 +« 4 ) 


and for Case (ii): 


<7i = />! +26 2 ; q 2 = b 3 ~bj 2 ; </ 3 = bj 3 


(5.7a) 

(5.7b) 

(5.8a) 


(5.8b) 


On the other hand, considering the definition of tensor M in Eq. (5. 2d), its functional dependence 
on C is of the same form as that of S in Eqs. (4.2) - (4.4) in terms of C for the unconstrained model. 
Thus, invoking the assignment “equations" (5.3), together with M(X^ <— S(Xj),, we can write 

(i) M = Eq. (4.2a) ; (ii) M = Eq. (4.2b) ; (iii) M = Eq. (4.2c) (5.9) 

where the principal values = M(\j) are given by (recall Eq. 5.1c) 

M(X,) = -l(A# fl +tf r t/ ) (5-10) 


Finally, by virtue of the continuity of the “new" tensors A, B, S, and M in terms of C, , and 
according to the previous analysis in Sec. 4, the expression for D in Eq. (5.2) is, ensured to be “well- 
defined" for the entire range of Xj. 

Remark 5.1 . In implementation, one actually uses C' 1 = (C 2 - Ij C + I 2 I)/l 3 - Also, explicit expressions 
are available for Xj and Xj directly in terms of C; e.g. using the (Ij J 2 , 0) set of invariants; see [17, p. 
269] for details. 
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6. COMPARISON WITH AVAILABLE RESULTS 


By way of a simple illustration, we compare some of the expressions obtained here with those 
available for the two-dimensional case in [e.g. 8]. To this end, we consider the case of (two-dimensional) 
double coalescence, and recollect the final observation (iii) made in Remark 4.2. Thus, for the model of 
Sec. 4, we simply use Eq. (4.8) here; i.e., with g = W ^ (same for any r), D = 4 g I^\ which, despite its 
extreme simplicity, leads to exactly the same coefficients; e.g., D nil = 4g, D 1122 = 0? and D 1212 = D 1221 
= 2g, etc., as those derived differently in [8,13]. 

Similarly, although not included here to limit the space, the detailed term-by-term expansions for 
D have also confirmed the equivalence of the forms given in [8,13], for the general two-dimensional case 
with distinct eigenvalues, to those obtained from Eq. (4.7) in the present three-dimensional case with 
double coalescence (with appropriate restrictions of the range of subscripts in S, C, E, etc.). Of course, 
modifications should be made in these comparisons to allow for different notations and definitions used 
(e.g. the engineering definition for shear strains in £ instead of the tensorial ones employed here, etc.). 

7. CONCLUSIONS 

Taken separately, the main constituents of the deformation tensor, i.e., principal values and 
associated eigenvectors, are, in general, not uniquely-defined and continuously differentiable functions. 
A careful consideration is thus called for in implementing constitutive models formulated in terms of these 
principal-strain measures, and this is the main problem addressed in this paper. In particular, the difficulty 
is entirely bypassed by resorting to explicit derivations of appropriate forms of the material tangent- 
stiffness matrices which are valid for the entire deformation range. These were developed here for two 
specific forms of the Ogden-type, strain-energy functions, which actually encompass many of the popular 
representations currently in use for rubber materials. A key feature in these is the underlying separability 
property, and this was employed to obtain the concise final forms of the tensor expressions, thus leading 
to their effective and robust numerical implementation. 

In view of the results obtained, it became obvious that the task simply amounts to application of 
a systematic limiting procedure for only one type of tensor-valued functions and their material time 
derivatives; i.e., those with symmetric single-argument functions in their spectral representation. In 
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particular, this was found essential in greatly simplifying the development given for the important case 
of uncoupled volumetric/deviatoric formulations. Finally, an important area for effective utilization of the 
present results would be to explore the feasibility of systematically deriving th e full expanded forms, 
together with generating the associated computer codes for implementation, using computer symbolic- 
manipulation packages [29]. 


8. REFERENCES 

1. A.E. Green and J.E. Adkins, Large Elastic Deformations, 2nd Edition, Oxford Univ. Press, 
Clarendon, 1972. 

2. R.S. Rivlin, “Large Elastic Deformations of Isotropic Materials; Part III: Experiments on the 
Deformation of Rubber", Philo. Trans. Roy. Soc . London, Ser. A , Vol. 243, pp. 251-288, 1951. 

3. R.S. Rivlin, “Forty Years of Nonlinear Continuum Mechanics", Proc. IX Int ., Congress on 
Rheology, Mexico, 1984. 

4. L.R.G. Treloar, The Physics of Rubber Elasticity , 3rd Ed., Oxford Univ. Press, U.K., 1975. 

5. R.W. Ogden, “Recent Advances in the Phenomenological Theory of Rubber Elasticity", J. Rubber 
Chem. Tech., Vol. 59, pp. 361-383, 1986. 

6. R.H. Finney and A. Kumar, “Development of Material Constants for Non-linear Finite Element 
Analysis", J. Rubber Chem. Tech .. Vol. 61, pp. 879-891, 1988. 

7. J.T. Oden, Finite Elements of Nonlinear Continua , McGraw-Hill, New York, 1972. 

8. T. Sussman and K.J. Bathe, “A Finite Element Formulation for Nonlinear Incompressible Elastic 
and Inelastic Analysis", Comp. Struct. , Vol. 26, pp. 357-409, 1987. 

9. J.C. Simo, R.L. Taylor, and K.S. Pister, “Variational and Projection Methods for the Volume 
Constraint in Finite Deformation Elastoplasticity", Comp. Meth. Appl. Mech. Engng.. Vol. 5 1 , pp. 
177-208, 1985. 

10. J.C. Simo and M. Ortiz, “A Unified Approach to Finite Deformation Elastoplasticity Based on the 
Use of Hyperelastic Constitutive Equations", Comp. Meth. Appl. Mech. Engng. , Vol. 49, pp. 
222-235, 1985. 


17 


11. S.T.T. Peng, “Nonlinear Multiaxial Finite Deformation Investigation of Solids’ Propellants", 
Report, AFRPI, TR-85-036, 1986. 

12. K.G. Valanis and R.F. Landel, “The Strain-Energy Function of a Hyper-elastic Material in Terms 
of Extension Ratios", J. Appl. Phys., Vol. 38, pp. 2997-3002, 1967. 

13. T.Y. Chang, A.F. Saleeb, and G. Li, “Large Strain Analysis of Rubber-Like Materials by a Mixed 
Finite Element”, Computational Mech., submitted, 1990. 

14. M.E. Gurtin, An Introduction to Continuum Mechanics. Academic Press, New York, 1981. 

15. J.E. Marsden and T.J.R. Hughes, Mathematical Foundations of Elasticity, Prentice-Hall, New 
Jersey, 1983. 

16. L.E. Malvern, Introduction to the Mechanics of a Continuous Medium, Prentice-Hall, New Jersey, 
1969. 

17. W.F. Chen and A.F. Saleeb, Constitutive Equations for Engineering Materials , Vol. 1 : Elasticity 
and Modeling, John Wiley & Sons, New York, 1982. 

18. T. Kato, Perturbation Theory for Linear Operators , 2nd Ed., Springer-Verlag, New York, 1976. 

19. P.R. Halmos, Finite-Dimensional Vector Spaces. 2nd Ed., Van Nostrand, Princeton, New Jersey, 
1958. 

20. J. Bathe, Finite Element Procedures in Engineering Analysis. Prentice Hall, Englewood Cliffs, 
New Jersey, 1982. 

21. Z.H. Guo, “Rates of Stretch Tensors", J. Elasticity , Vol. 14, pp. 263-267, 1984. 

22. M.E. Gurtin and K. Spear, “On the Relationship Between the Logarithmic Strain Rate and 
Stretching Tensor", Int. J. Solids Struct. , Vol. 19, pp. 437-444, 1983. 

23. A. Hoger and D.E. Carlson, “On the Derivative of the Square Root of a Tensor and Guo’s Rate 
Theorems", J. Elasticity, Vol. 14, pp. 329-336, 1984. 

24. R. Hill, “Constitutive Inequalities for Isotropic Elastic Solids under Finite Strain", Proc. Roy. Soc. 
London, Vol. A3 14, pp. 457-472, 1970. 

25. R. Hill, “On Constitutive Inequalities for Simple Materials - Part I", J. Mech. Phys. Solids, Vol. 
16, pp. 229-242, 1968. 


18 


26. J.W. Hutchinson and K.W. Neale, “Finite Strain J2 Deformation Theory", in Finite Elasticity, D.E. 
Carlson and R.T. Shield (Eds.), Martinus Nijhoff Publisher, The Hague, The Netherlands, pp. 
237-247, 1982. 

27. W.D. Rolph III and K.J. Bathe, “On a Large Strain Finite Element Formulation for Elastoplastic 
Analysis", in Constitutive Equations: Macro and Computational Aspects. K.J. Wiliam (Ed.), 
American Soc. Mech. Engng., New York, pp. 131-147, 1984. 

28. J.M. Ball, “Differentiability Properties of Symmetric and Isotropic Functions". Duke Mathematical 
Journal. Vol. 51, pp. 699-728, 1984. 

29. S.M. Arnold and H.Q. Tan, “Symbolic Derivation of Potential Based Constitutive Equations", 
Computational Mech., Vol. 6, pp. 237 - 246, 1990. 

30. J. Betten, “Applications of Tensor Functions to the Formulation of Constitutive Equations, 
Involving Damage and Initial Anisotropy", J. of Eng. Fracture Mech.. Vol. 25, No. 5/6, pp. 573- 
584, 1986. 


19 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for re 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments reg< 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget. Paperwork Reduction 1 

^viewing instructions, searching existing data sources, 
irding this burden estimate or any other aspect of this 
information Operations and Reports, 1215 Jefferson 
Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY ( Leave blank) 

2. REPORT DATE 

November 1991 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 



5. FUNDING NUMBERS 


REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Explicit Robust Schemes for Implementation of a Class of Principal Value- 
Based Constitutive Models: Theoretical Development 


6. AUTHOR(S) 

A.F. Saleeb and S.M. Arnold 


WU-510-01-50 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-6714 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 105 345 


11. SUPPLEMENTARY NOTES 

A.F. Saleeb, University of Akron, Department of Civil Engineering, Akron, Ohio 44325; S.M. Arnold, 
NASA Lewis Research Center. Responsible person, S.M. Arnold, (216) 433-3334. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 39 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

The issue of developing effective and robust schemes to implement a class of the Ogden-type hyperelastic constitu- 
tive models is addressed. To this end, explicit forms for the corresponding material tangent stiffness tensors are 
developed, and these are valid for the entire deformation range; i.e., with both distinct as well as repeated principal- 
stretch values. Throughout the analysis the various implications of the underlying property of separability of the 
strain-energy functions are exploited, thus leading to compact final forms of the tensor expressions. In particular, this 
facilitated the treatment of complex cases of uncoupled volumetric/deviatoric formulations for incompressible 
materials. The forms derived are also amenable for use with symbolic-manipulation packages for systematic code 
generation. 


14. SUBJECT TERMS 

Constitutive equations; Isotropic; Eigenvalves; Eigenvectors; Stiffness matrix 


15. NUMBER OF PAGES 

20 


16 . 


PRICE CODE 

A03 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 


